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ABSTRACT 

We present a systematic study of the steep decay emission from gamma-ray 
bursts (GRBs) observed by the Swift X-Ray Telescope (XRT). In contrast to 
the analysis described in recent literature, we produce composite Burst Alert 
Telescope (BAT) and XRT light curves by extrapolating the XRT data (2-10 
keV) into the BAT energy range (15-25 keV) rather than extrapolating the BAT 
data into the XRT energy band (0.3-10 keV). Based on the fits to the composite 
light curves, we have confirmed the existence of an exponential decay component 
which smoothly connects the BAT prompt data to the XRT steep decay for 
several GRBs. We also find that the XRT steep decay for some of the bursts can 
be well fit by a combination of a power-law with an exponential decay model. We 
discuss this exponential component within the frame work of both the internal 
and the external shock model. 
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Introduction 



The transition between the gamma-ray burst (GRB) prompt emission and afterglow 
emission has generated great interest within the scientific community. It is generally ac- 
cepted that the GRB prompt emission is due to internal shocks originating from the col- 
lision of faster and slower moving shells, whereas, the afterglow is believed to originate 
from an external shock resulting from the relativistic fireball colliding; with a circum-burst 



medium ( jRees fc Meszarod Il994 ; ISari fc Piranl 119971 ; iMeszaros fc Reed 119971 ) . During the 



GRB episode there should be a transition from one phase to the other, however, it is still 
not well understood as to when this transition occurs. The Burst and Transient Source 
Experiment (BATSE) observation of GRB 980923 showed l ong lasting tail em ission, ~400 



s, which is best described by a power-law temporal de cay (jGiblin et al.lll999l) . Based on 
the spectral and temporal characteristics of this burst, Giblin et al. ( 19991 ) concluded that 



the tail emission was a part of the after glow emission, thus , the external shock could be 



generated during the prompt 7-ray phase. iGiblin et al.l (120021 ) performed a systematic study 



of the prompt tail emission using 40 GRBs observed by BATSE. They found that the tem- 
poral decays are best described by a power-law with a decay index of —2 rather than an 
exponenti al. There are several ot her analysis of BATSE data which reach the same conclu- 



sion (e.g.. lRyde fc Svenssonll2002l ). According to the BeppoSAX observations, the late time 



afterglow smoothly connects with the prompt emission if the onset time of the light curve 



Pian et al. 


2001; 


Piro et al. 


2005) 



(e.g., 

X-ray spike represents the onset of the afterglow. However, the delay of a few hours to a few 
days before the narrow field X-ray instrument is pointed to the GRB position, weakens the 
discussion concerning the transition from the prompt emission to the afterglow. 



As a result of the revolutionary features of Swift (IGehrels et al.ll2004l ). our understanding 
of the X-ray properties of GRBs has been improved dramatically. With the combination of 
the accur ate on-board calculation of the GRB position by the Burst Alert Telescope (BAT: 15- 
150 keV; iBarthelmy et al.l 120051 ) and the fast slewing capability of the spacecraft, Swift 
can begin a highly s ensitive X-ray observation with the X-Ray Telescope (XRT:0.2-10 keV; 
Burrows et al.ll2005l ) within a few tens of seconds to a few hundred seconds after the burst 
trigger. According to t he XRT observations, the X-ray proper ties of the GRB emission have 



very complex features (INousek et al.ll2006l ; IZhang et al. 



2006al ). One of the most unexpected 



discoveries by the XRT is the existence of the steep decay component during the initial phase 
of the X-ray light curve. The origin of this steep decay component is generally considered to 
be a result of the delayed pro mpt emission from different viewing latitudes of the jet, the so 



called "c urvature effect" (e.g. iFenimore et al.lll996l; iKumar fc Panaitescull2000l ; IZhang et al. 



2006a ). Tagliaferri et al. ( 2005 ) and Barthelmy et al. ( 20061 ) investigated the steep decay 
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component with the composite BAT and XRT light curves for several GRBs. To generate 
the composite light curve, both papers performed an extrapolation of the BAT mask- weighted 
(background subtracted) light curve into the XRT 0.2-10 keV energy band using a best-fit 
power-law photon index. The authors found that GRB 050126 and GRB 050219A do not 
show continuous emission from the BAT to the XRT light curve, however, GRB 050315 
and GRB 050319 do display a smooth continuation from the BAT to the XRT light curve. 
O'Brien et al.l (120061 ) performed a systematic study of the early X-ray emission using a sample 
of 40 Swift GRBs. They constructed a composite BAT and XRT light curve in the 0.3-10 
keV band. In order to extrapolate the BAT data points, the BAT mask-weighted count rate 
was converted to flux in the 0.3-10 keV band using the mean of the best-fit photon indices 
obtained from a simple power-law fit to both the BAT and the XRT spectrum. The authors 
found that a fit of the 40 superimposed GRB light curves in the 0. 3-10 keV band could be 
described by an exponential decay followed by a power-law decay. IWillingale et al.l (120071 ) 
investigated the X-ray emission including the late time light curve data. They found that 
the X-ray light curve data can be modeled with the superposition of an early ("prompt") 
and a late time ("afterglow") component. 

In this paper, we describe the analysis of the tail emission using an alternative ap- 
proach. We generate the composite BAT and XRT light curve in the 15-25 keV BAT energy 
band by extrapolating the 2-10 keV XRT count rate into the 15-25 keV band. Either 
approach, extrapolating the BAT count rate down to the XRT band (hereafter BAT-to- 
XRT extrapolation) or extrapolating the XRT count rate up to the BAT band (hereafter 
XRT-to-BAT extrapolation), can encounter similar syst ematic problems. Most of the GRB 
prompt emission spectra is well fit by the Band function (IBand et al.lll993l) with a low-energy 
and a high-energy photon index of ~ 1 and ~ 2.3, respectively (e.g. ISakamoto et al.l 12005 



Kaneko et al.ll2006l ). However, it is a well known chara cteristic that the peak energy of the 



spectrum, E pp „ k , not only evolves during the burst fe.g.. lLlovd-Ronning and Petrosian 112002 



Frontera et al.lll999l ) but also changes from burst-to-burst (e.g.. lSakamoto et al.ll2005l ). Here, 



we discuss the issues with regard to both types of extrapolation by considering five cases 
depending on the value of E peak (shown in figure [1] left is the BAT-to-XRT extrapolation 
and right is the XRT-to-BAT extrapolation). If we consider the energy bands of the BAT 
(15-150 keV) and the XRT (0.3-10 keV), the five cases are defined by: 1) E pcak > 150 keV, 
2) 15 keV < £ pcak < 150 keV, 3) 10 keV < E pcak < 15 keV, 4) 0.3 keV < E peak < 10 keV, 
and 5) E peBk < 0.3 keV. For cases 1 and 5, both the BAT-to-XRT and the XRT-to-BAT 
extrapolations should provide the correct flux because the photon index of the extrapolated 
energy band is the same as the observed energy band ("1" and "5" of the left and the right 
panels of figured]). For cases 3 and 4, the BAT-to-XRT extrapolation would over estimate 
the flux in the XRT energy band, since the photon index in the observed energy band is 



-4- 



steeper than that of the extrapolated energy band ("3" and "4" of the left panel of figure 
[1]). Similarly, for case 2 and 3, the XRT-to-BAT extrapolation would over estimate the flux 
in the BAT energy band, because the photon index in the observed energy band is shallower 
than that of the extrapolated energy band ("2" and "3" of the right panel of figured]). For 
case 2, using the BAT-to-XRT extrapolation ("2" of the left panel of figured]), and for case 
4, using the XRT-to-BAT extrapolation ("4" of the right panel of figured]), the flux would 
be overestimated in the XRT and the BAT energy band respectively, because the simple 
power-law fit is a tangential line to the curved spectrum as a result of the narrow energy 
bands of the BAT and the XRT (e.g. Sakamoto et al. in preparation). 

To date, most of the BAT and XRT composite light curves in the literature have been 
produced from a BAT-to-XRT extrapolation. However, we believe that the XRT-to-BAT 
extrapolation described here may minimize the systematic errors, especially with respect to 
investigating GRB tail emission, for the following three reasons. First, we can account for the 
spectral evolution during the prompt emission if we extract the flux from the time-resolved 
spectral analysis of the BAT data. As previously mentioned, -E pea k shifts from hundreds of 
keV to a few keV during the prompt emission. If in the analysis, one does not account for 
the spectral evolution during the prompt emission, which is the case for the majority of the 
published XRT and BAT composite light curves using the BAT-to-XRT extrapolation, the 
systematic error in the extrapolated flux in the X-ray range could be significant. Second, 
since the BAT mask- weighted count rate does not correct for the energy dependence of each 
photon, the count rate of the source in the off-axis direction will be systematically smaller 
than the on-axis case. This effect becomes an issue when the BAT-to-XRT extrapolation 
has been performed by converting the BAT mask-weighted count-rate into flux using a fixed 
photon index obtained from a simple power-law fit. According to the BAT Crab observation, 
the count rate from the Crab is ~15% smaller in the 45 degree off-axis case. This off-axis 
effect is correctly taken into account in the BAT energy response matrix, but not in the BAT 
mask-weighted count rates. Therefore, unless one applies an additional off-axis correction to 
the BAT mask-weighted count rates, a systematically smaller flux is obtained if the source 
is in the off-axis direction, which is always the case for the BAT GRB data prior to the 
spacecr aft slewing to t he GRB position. Third, according to the GRB synchrotron shock 



model (jSari et al.lll998l ). if the observed spectrum has a photon index steeper than 2, and a 
power-law index of an electron distribution, p, in the range of 2 < p < 3 (where N(7 e ) oc % p , 
and 7 e is the Lorentz factor of the electrons), then the observed frequency should be above 
the synchrotron critical frequency for electrons with a minimum Lorentz factor (u m ) in the 
fast cooling phase or above the cooling frequency (z/ c ) in the slow cooling phase. In this case, 
since there is no characteristic frequency above u m (in the fast cooling case) and v c (in the 
slow cooling case), it is reasonable to extrapolate upwards in energy. The electron power-law 
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index 2 < p < 3 is typic al for both the prompt em i ssion (e.g., 
afterglow of GRBs (e.g.. IPanaitescu k, Kumarj|2002l ; lYost et al. 



Kaneko et al.ll2006l ) and the 



20031 ). Thus, if we select the 



burst samples which have an observed photon index steeper than 2 in the XRT data, the 
systematic error in the XRT-to-BAT extrapolation should be reduced. 

Following these arguments, we describe in this paper the XRT-to-BAT extrapolation to 
investigate the nature of the XRT steep decay component. As we mentioned, in principle the 
XRT-to-BAT extrapolation has similar issues to the BAT-to-XRT extrapolation. However, 
for certain bursts the XRT-to-BAT extrapolation could greatly reduce the systematic error. 
More importantly, a different approach to the analysis may provide an alternate view of the 
problem. 



2. Analysis 

X-ray light curves were produced from the XRT data for all Swift GRBs detected be- 
tween June 2005 and September 2006. 54 GRBs with an early phase power-law decay 
index steeper than —2 (to taken as the BAT trigger time) were identified. From the spectral 
analysis of the X-ray data from these bursts in the 0.3 - 10 keV band, we found that many 
bursts exhibited a significant difference in the best-fit spectral parameters between the early 
(< 1000 seconds) data and the later (> 1000 seconds) data. This is likely due to the fact that 
-E"peak is moving through XRT energy range early in the observation but one cannot measure 



this f rom the XRT data alone because of the narrow energy band (e.g., iButler fc Kocevski 



20071 ). To minimize the effects of the spectral evolution and also the absorption, only data 
above 2 keV was used in these analysis. We further refined the burst sample to those which 
satisfied the following criteria: 1) Greater than four data bins containing at least 20 counts 
in the 2-10 keV light curve during the early phase (< 1000 seconds). This criteria reduced 
the number of bursts in the sample significantly, due to much fewer counts being detected 
above 2 keV. 2) No more than a single flare present in the early XRT light curvqj. 3) The 
joint spectral analysis of the Photon Counting (PC) and Windowed Timing (WT) data must 
include a photon index of 2 in the 90% confidence interval. 13 GRBs satisfied the screening 
criteria. Since the 2-10 keV joint spectral fit to the PC and WT data includes a best fit 
photon index steeper than 2 for all the GRBs in the sample, -E pea k should be below 2 keV. 



1 Following updates to the on-boar d software which com pensate for the uncontrolled temperature due to 



failure of the co oling control system (jKennea et al. 

(jam et aiJIiooih . 

2 Bursts exhibiting flares in the BAT were not excluded. 



20051 ). bright Earth effects and micrometeriod damage 
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This will reduce the systematic error of the extrapolation because there should be no spectral 
evolution as a result of the shift of E pe3ik during the burst. 



As mentioned in IZhang et al.l (l2006al ). the definition of the offset time (to) is critical 
when performing fits to the early phase light curve. Traditionally, to is defined as the trigger 
time of the GRB instrument; when the count rate exceeds some background level (rate 
trigger). However, the definition of the BAT trigger time is different. The BAT trigger time 
is the start time of the foreground time interval of the image from which the GRB is detected 
on-board. Thus, to be comparable to a rate trigger time, we define to as the start time of 
the prompt emission (start time of tioo interval) for the whole sample. 



2.1. BAT analysis 

The BAT analysis was performed using HEAsoft version 6.1.1 and CALDB version 
2006-05-30. The event-by-event data were used for these analysis. The non- linear energy 
correction for each event was applied by bateconvert. The mask- weighting factors were 
calculated by batmaskwtevt using the on-board position. The detector enable and disable 
map was created by bathotpix combining the enable and disable map generated by the 
flight software. We created the BAT light curve by batbinevt in the full energy range 
(15-350 ke\fl) in 4 ms bin except for GRB 051109A (64 ms), GRB 060427 (1 s) and GRB 
060923C (64 ms). We used larger binning for these three GRBs because of the low signal- 
to-noise ratio of the emission. The duration and the time intervals based on the Bayesian 



Block algorithm (jScargld 119 98) were calculated by battblocks. The spectrum of each time 
interval was extracted by batbinevt. The energy response file was created by batdrmgen. 
If the time interval was during the spacecraft slew, we updated the keywords in the spectral 
file related to the energy response process by batupdatephakw and then created the energy 
response file for the time interval by batdrmgen. We applied systematic error vectors to the 
spectrum using batphasyserr prior to doing the spectral analysis. The spectral analysis 
was performed using Xspec 11.3.2. 

The energy flux in the 15-25 keV band was calculated for each time interval directly from 
the spectral fitting process. The spectra from each time interval were fitted with a simple 
power-law model. According to the BAT GRB catalog (Angelini et al. in preparation), the 
detection threshold of the BAT in the 15-25 keV band is ~ 10~ 9 ergs cm -2 s _1 . Based 
on this result, the 15-25 keV flux was treated as an upper limit when the calculated 15- 



3 The coded mask is transparent to photons above 150 keV. Thus, photons above 150 keV are treated as 
background in the mask-weighted method. The effective upper boundary is ~ 150 keV. 
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25 keV flux was less then 10 -9 ergs cm -2 s~ x . The upper limit was estimated from using 
the event-by-event data from the Crab nebula on-axis observation collected on 2005 March 
24 (observation ID: 00050100016). According to this observation, the BAT can detect the 
Crab nebula in the 15-25 keV band at 5 a in a one second e xposure. Assuming tha t the 



BAT sensitivity scales as the square-root of the exposure time ( iMarkwardt et al.l 120051 ) and 
a canonical Crab flux of 5.3 x 10 -9 ergs cm -2 s" 1 in the 15-25 keV band, we calculated the 
BAT 5 a upper limit in the 15-25 keV band from the following equation, 

K Q x 1 n -9 

F(15 - 25keV) 5CT = tj^ (ergscm" 2 s" 1 ). (1) 

Jpcode 

Here t exp is the exposure time and f pco de is the partial coded fraction. Since our estimation 
of the 5<t upper limit is based on the on-axis Crab observation, f pco de will correct for a source 
observed in the off-axis direction. 

For the time-averaged spectral analysis, we use the time interval from the emission start 
time to the emission end time (t±oo interval). When the spacecraft slew occurred during 
the time interval, we created the response matrices for each five second period, taking into 
account the position of the GRB in detector coordinates. We then weighted these response 
matrices by the five second count rates and created the averaged response matrices. Since 
the spacecraft slews about one degree per one second in response to a GRB trigger, we choose 
five second intervals to calculate the energy response for every five degrees. 



2.2. XRT analysis 



The 13 bursts meeting the criteria outlined in section 2 were processed using the HEAsoft 
tools version 6.1.1 including version 2.5a of the Swift software. The level 2 cleaned event 
files were produced from the xrtpipeline task (version 10.4) using the standard screening 
criteria. Version 8 of the response matrices in CALDB and the co rresponding ancillary 
response produced from xrtmkarf were used. The standard grades (IBurrows et al.l 120051 ) 
were used in the analysis; grades 0-2 and 0-12 for WT and PC mode, respectively. 



Both WT and PC mode data (jHill et al.ll2004l . 120051 ) from the first observation segment 
(000) were analyzed. The source and background extraction regions were nominally a 40- 
pixel square and an annulus with a 3 pixel inner radius and 25 pixel outer radius for WT 
and for PC mode, respectively. The extraction region s were modified for the piled-up cases 
in WT mode in accordance with lRomano et al.l (120061 ). eliminating an inner square centered 
on the source. It is well documente d that pile-up in PC mode causes a redistribution of 
single pixel events to higher grades (jVaughan et al.ll2006l ). In order to account for pile-up 
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in PC mode, the grade distributions were analyzed for each burst during periods where the 
count-rate exceeded 0.5 counts/sec. For each burst, the percentage of single pixel events were 
plotted versus the radius of the center of the annulus for increasing radii (0-8 pixels). The 
annulus with smallest radius at which the percentage of single pixel events became constant 
was used for the PC mode analysis. Nominally, for a spectrum where the average energy is 
1.5 keV, 78% of the events will be single pixel events (IMoretti et al.l 12004 ) . 

The hardness ratio ((2. 0-10. 0)keV/(0. 5-2.0) keV) was examined for each burst to en- 
sure that there was no significant spectral evolution. An exposure map was created from 
xrtexpomap to correct for the dead columns and hot pixels, xrtmkarf was used to create 
the ancillary response files. This includes corrections for losses in the wings of the point 
spread function and the center of the annulus, for the exposure and for vignetting. 

The spectrum file was binned with a minimum of 20 counts/bin in order for x 2 statistics 
to be valid for the spectral fitting. Xspec version 11.3.2 was used to perform a joint spectral 
fit of the WT and PC data from 2-10 keV using a simple power law model. The 15-25 keV 
normalization obtained from the pegpwrlw model was used to extrapolate the XRT count 
rate into flux in the BAT 15-25 keV energy range. Only light curve bins with greater than 
90% exposure were used in order to limit errors due to dead-time incurred by instrument 
mode switching. The light curve was binned to have greater than 20 counts/bin. 



2.3. Fitting a composite BAT and XRT light curve 

To investigate the connection between the prompt and the afterglow emission in the 
composite light curves, we first fit the XRT light curve, only using data where the hardness 
ratio showed no strong spectral evolution. We then fit both the BAT and XRT light curves 
jointly. Both fits were performed using a power-law model with an offset time (PLO), 

^15-25keV = K pow (t — t P Q W ) a , (2) 

where t^ w is a offset time, a is a decay index, and K pow is the normalization, and with an 
exponential model (EXP), 

^i5-25kcv = K exp exp( ), (3) 

w 

where w is the decay constant and K exp is the normalization. For the XRT only fit, we 
fixed to° w to zero. Finally, we fit the BAT and XRT light curves simultaneously using a 
combination of a power-law model with an exponential decay component (PLEXP), 

i 71 i5-25kev = K pow (t - tl° w )~ a + K exp exp( ). (4) 
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For the fit to the XRT only light curve, the time interval was from the first XRT data 
point to the last data point before showing a deficit from the PLO model using an offset 
time of zero. For the joint BAT and XRT fit, the time interval was from the first BAT 
data point to the last XRT data point before showing a residual from the PLO model. Any 
other definitions used for the time interval are stated as a footnote in table [D The values of 
the time intervals used in the light curve fitting are shown in the fourth column of table [TJ 
The best-fit model was selected based on the x 2 °f the fit. However, because a PLO model 
will not fit the data before t^ w , the judgment between PLO and PLEXP is based on visual 
inspection as to whether the model fit both the BAT and XRT data simultaneously or not. 



3. Results 

The left panels of figures [2H6] show the composite BAT (black open circles) and XRT 
(red open triangles) light curve in the 15-25 keV band overlaid with the best-fit light curve 
model. The light curve models are PLO (eq. ©), EXP (eq. ©), and PLEXP (eq. (g])) 
from top to bottom. In the bottom figure of PLEXP, both PLO and EXP components are 
also shown as a dashed and a dash-dotted line, respectively. The best-fit parameters of the 
light curves are summarized in table [D The best-fit of these models is labeled in blue. The 
right panels of figures [2HH] from top to bottom show the BAT light curve in the 15-150 keV 
band, the BAT photon index based on the time-resolved spectral analysis, the XRT 2-10 
keV count rate, and the XRT count rate ratio (2.0-10.0) keV/(0.5-2.0) keV. The best-fit 
spectral parameters based on the 2-10 keV joint fit to the XRT WT and PC mode data are 
summarized in table [2j 

From the initial steep decay phase of the XRT light curve, it is difficult to distinguish 
between PLO and EXP from the XRT data alone. Both models fit equally well for all of the 
bursts. However, the difference and the importance of the individual components become 
clear when the BAT data are included in the fit. An EXP model fits well for GRB 050814, 
GRB 0509 15B, GRB 060427 and GRB 060428B. A PLO or a PLEXP model is not required 
for these GRBs. For GRB 060923C, a PLO is the model best represented by the composite 
light curve. A PLEXP is the best model for the remaining 7 GRBs. The best-fit parameters 
which we used in the systematic study presented in this section, are shown in bold font in 
table dJ 



First, we investigated the possibility of the curvature effect for those GRBs which have a 
PLO c omponent in the composite BAT and XRT light curve fit. According to the curvature 



effect (IFenimore et al.l Il996l ; iKumar &: Panaitescul l2000i ; IZhang et al.l l2006al ). the relation 
between the decay index, a, and the XRT photon index, T X rt, should be described by, a 
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= 1+Txrt, if the curvature effect is the cause of the XRT steep decay. Figure [7] shows the 
correlation between a and Txrt for our sample. The dashed line is the expected relationship 
from the curvature effect (a = 1+Y X rt)- Although GRB 060923C may be consistent with 
the curvature effect, the majority of the bursts in the sample do not satisfy the expected 
relation. The inconsistency with the curvature effect could be due to neglecting the spectral 
evolution during the steep decay in our analysis of the XRT data. Looking at the time 
evolution of the count rate ratio of our sample, we find a hard to soft evolution from 0.6 to 
0.5 and from 0.6 to 0.4 for GRB 060202 and GRB 060211A. These changes correspond to an 
evolution of the photon index from 1.5 to 1.7 and from 1.5 to 1.9, respectively, according to 
the calculation by the Xspec fakeit command using the detector and the ancillary response 
files created for each source region. If this spectral evolution is taken into account, the steep 
decay could be consistent with the expectation of the curvature effect for GRB 06021 1A and 
GRB 060202. However, we do not see a strong spectral evolution for the other GRBs with 
the exception of GRB 060418. Note that for GRB 060418, there is a strong evolution in the 
hardness ratio during the episode at to+150 s which may cause an error in the extrapolated 
flux. 

We find that most of the sample requires an EXP component to fit the BAT and XRT 
light curves simultaneously. Therefore, we can conclude that some of the early steep decay 
observed by XRT is a continuation of the exponential decay tail of the prompt emission. 
Interesting characteristics can be found for the bursts where a PLEXP model is the most 
representative model for the composite light curve. The dominant component of the fit to 
the XRT light curve 180 s after t for GRB 060202 is an EXP. For GRB 050803 and GRB 
060109 there is almost equal contribution from the EXP and PLO components in the initial 
XRT data . Whereas, a PLO is the dominant component for GRB 051109A, GRB 060111B 
GRB 06021 1A, GRB 060306 and GRB 060418B. This result clearly demonstrates that the 
XRT steep decay could be composed of at least two different components. Without careful 
consideration of both the BAT and the XRT data simultaneously, it is not p ossible to dis 



tingui sh between these two different components. It is important to note that lO'Brien et al. 



( 120061 ) also reached a similar conclusion; that the BAT and XRT composite light curve is 



composed of an exponential decay which relaxes to a power-law decay. 

For the bursts that exhibit an EXP component, we investigated the correlation between 
the exponential decay index, w, and the prompt emission properties derived from the BAT 
data (table [3]). The results are summarized in figure No correlation is found for the 
properties of the prompt emission except between w and the BAT T 90 which is expected 
because both parameters are related to the duration of the bursts. 

Based on our study, there is a strong indication that the steep decay component observed 



- 11 - 



by the XRT is part of the prompt emission (e.g., also lNousek et al.ll2006l ; lQ'Brien et al.ll2006l ). 
Thus, we calculate the fluence which may be below the sensitivity limit of the BAT. This 
fluence was calculated by accumulating the flux from the best-fit composite BAT and XRT 
light curve model from the end of the emission as detected by the BAT to 1000 seconds after 
to- Figure [9] shows the ratio of the percentage of the fluence in the tail emission which is below 
the BAT sensitivity limit and the fluence recorded by the BAT. For 7 out of the 13 GRBs 
in the sample, the fluence of the tail component is less than 15% of the fluence recorded by 
BAT. However, more than 15% of the fluence may be radiated below the BAT sensitivity 
limit for GRB 050915B, GRB 051109A, GRB 060202, GRB 060211A, GRB 060427, and 
GRB 060428B. This result gives rise to the question as to whether the fluence measured by 
the 7-ray instrument reflects the true fluence of the prompt emission. 



4. Discussion 

We have presented the BAT and XRT composite light curves, derived extrapolating the 
XRT 2-10 keV flux up to the BAT 15-25 keV energy range for GRBs which have a steep 
decay component in the initial XRT light curve. Based on the simultaneous fit of both the 
BAT and XRT light curves, we have confirmed the existence of an EXP component which 
smoothly connects the BAT prompt emission to the XRT steep decay for several GRBs. 
We have also found that the XRT steep decay for some of the bursts can be fit well by a 
PLEXP model. In the following sections, we discuss the possible origins of the PLO and 
EXP components. 



4.1. Origin of the PLO component 



A PLO component most likely originates from an internal shock (so called, curvature 
radiation or high- latitude emission associated with the last bright sp ike. iKumar fc Panaitescu 
2000l ; iNousek et al.ll2006l ; IZhang et al.ll2006at lYamazaki et al.ll2006l ) . Our results support this 
idea, because for most bursts the XRT steep decay component smoothly connects with the 
last bright episode detected by the BAT (e.g. GRB 050 803). The instantaneous em ission 
from a uniform jet produces a decay index of a = 1 + T fcumar fc Panaitescul lioool ). This 
formula was examined using the power-law decay index derived from the PLO model and 
the photon index based on the joint WT and PC spectral analysis of our sample. For the 
majority of the sample, we find that the decay index is not consistent with the formula. One 
of the possible reasons for this inc onsistency could be the spectra l evolution during the early 
XRT observation in some bursts (IZhang. Liang, fc Zhang||2006l ). However, as discussed in 
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section 3, the inconsistency cannot always be associated with spec t ral evolution. Anothe r 
possible rea s on is the choice of the time-zero (jZhang et al.l l2006al ; lYamazaki et al.l 120061 ) . 



Liang et al.l (120061 ) i nvestigated the curva ture effect as an origin of the XRT steep decay 



using the data set of lO'Brien et al.l (120061 ). They made the assumption that the XRT steep 
decay component is due to the curvature effect and investigated whether the time-zero is 
consistent with the beginning of the bright episode. They concluded that for most of the 
sample, the time-zero was consistent with this picture. The main issue with their approach 
is that for the fixed power-law decay index, as expected from the formula a = 1 + T, by 
shifting the offset time it is possible to fit the decay index to most early XRT data. This is 
because in their approach, the fitting parameters are not only the offset time but also the 
normalization, allowing extra freedom in the fit. Here, we demonstrate this problem using 
GRB 050803 and GRB 050814. The dash-dotted lines in figure M show the best fit PLO 
model by changing a from 1 to 5 for GRB 050803 and from 2 to 3 for GRB 050814, and 
va rying tn° w . Only th e XRT data, shown in red triangles, are used in the fitting process as 



m 



Liang et al.l (120061 ). As seen in the figures, the choice of a and t^"" is not unique if one 



only try to fit the XRT data. Moreover, as it is clearly demonstrated in the case of GRB 
050814, even if t^"" is chosen as the start time of the GRB pulse, the intensity of the pulse 
expected from the model is an order of magnitude brighter than the data. Therefore, if a 
bright episode in the BAT data, which could contribute to the steep decay component, is 
simultaneously fit with the XRT data, as in our approach, both the offset time and the decay 
index will be uniquely constrained by the data. It may be difficult to test the curvature effect 
or the relation a = 1 + T definitely without fitting the XRT and BAT data simultaneously. 
The third possibility is to abandon the assumption of the uniform jet emission. Our results 
may suggest that the structure of the jet is much more complex than a uniform jet (e.g. 



Yamazaki et al.ll2006l ). 



4.2. Origin of EXP component 

Since it is difficult to explicitly state the origin of an EXP component, we will discuss 
the possibilities of both an internal shock and an external shock as the origin of an EXP 
component. 



4-2.1. External shock scenario for EXP component 



One interpretation of an E XP component is the prese n ce of the external shock emis- 
sion during the prompt phase. iFenimore &: Ramirez-Ruizl (119991 ) studied the case of the 
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co-existence of the emission from the external shock (deceleration of the initial shell) during 
the emission from the internal shocks. They showed that the smooth long lasting soft emis- 
sion which ar ose from an external shock could o verlay the light curve of the prompt emission. 
Furthermore, iFenimore &: Ramirez-Ruia (119991 ) showed that the efficiency of converting the 
bulk energy to radiation is 85% in this case, whereas, internal sh ocks without deceler ation 
only convert about 1% (IKumar I Il999c iPanaitescu et al.lll999l ). IZhang et al.l (j2006bl ) cal- 
culated the radiation efficiency using the Swift X-ray afterglow data and show that about 
half of the GRBs have an efficiency > 1% which provides a challenge for producing the 
prompt emission from internal shocks alone. Furthermore, according to the optical observa- 
tion by Rapid Telescopes for Optical Response (RAPTOR) during the prompt emission of 
GRB 05082 0A, smoothly decayin g emission which does not correlate with the prompt spikes 
was found (jVestrand et al.ll2006l ). If this emission is from the external shock, we might be 
observing the deceleration of the outflow during the prompt phase. 



Following the argument of IFenimore et aJj (119961 ) and lFenimore fc Ramirez- Ruizl (119991 ). 
we estimate the bulk Lorentz factor of our sample, assuming an EXP component is purely due 
to external shock emission. Let us assume that an external shock starts its emission at the 
radius, -R - The external shock will be decelerated because of sweeping up the inter-stellar 
medium (ISM). The total energy of the central engine, Eq, can be expressed as, 



E = (47r/3)i?onisMm p c 2 7o 



(5) 



where tt-ism is a density of the ISM, m p is the proton mass, c is the speed of light, and 70 is 
the bulk Lo rentz factor at i?o- The duration of the emission (AT) is determined by the radial 
time scale (jPiranl |l999_|) which is the difference in the arrival time of the photons emitted 
between R and ai? (a > 1) as measured by the observer, 

AT = [(a 4 -l)/4](i?o/2 7o 2 c). 

Therefore, the bulk Lorentz factor can be expressed as, 

70 = (327rm p c 2 /3)^/ 8 [(a 4 - l)/Af/\E,/n ISM ) l l\c^T)-^. (6) 

The relationship between AT and the full width at half maximum (FWHM) of the pulse 
(T1/2) is described by, 0.22 AT = T/ 2 , which is valid for a pulse shape of a fast rise, 



exponential decay (FRED) (IFenimore et al.l 119961 ). Thus, once we know the redshift, since 
T1/2 can be estimated from the best-fit parameters of an EXP component, we can calculate 
the bulk Lorentz factor as a function of E /tiism- I n the following arguments, we use a 
typical value of 2 for the parameter "a" (in the a = 2 case, AT is the radial time scale from 
R to 2Rq). 
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For bursts with unknown redshifts (all bursts in the sample except GRB 050814, GRB 
051109A and GRB 060418), we used the mean redshift of 2.4 obtained from the Swift long 



We can derive a reasonable range of 70 from 143 to 350 and Rq of ~ 10 16 cm 



GRBsfl. 

assuming .E 0/ Mism = 1 x 10 52 erg cm 3 for the 12 GRBs in our sample which have an EXP 
component in the best fit model (table [3]). According to the calculation of 70 by the Rapid 
Eye Mount (REM) telescope using th e peak time of the ear ly afterglow data, 70 is about 400 
for GRB 060418 and GRB 060607A jMolinari et al.ll2006h . Their 70 value agrees within a 
factor of two of our estimates. For GRB 060418, E Q will be ~ 1 x 10 54 ergs if we assume 70 
of 400 as derived from the REM observatio n and ?t.tsm of lcm" 3 . Eq of ~ 10 54 ergs is also a 
typical value according to the calculation of IZhang et al.l (j2006bl ) using X-ray afterglow data 
observed by XRT. 

It would be interesting to look for a correlation between 70 and E pea ^ in the GRB rest 
frame (i^ak) although it is difficult to calculate for our sample because we do not have 
measurements of both -Ep ea k and redshift. One of the advantages of using the Swift sample 
for this study is that soft GRBs, so called X-ray Flashes (XRFs), are included in the sample 
because of the relati vely softer energy r esponse of BAT compared to that of BATSE. In the 
dirty fireball model (IDermer et al.lll999l ). E^ k has a strong dependency on the bulk Lorentz 
factor (Ep f y oc 7 4 ). The un ified jet models for XRFs and GRBs, such as the structured- 



jet model (|Rossi et al.ll2002l ). and the variable jet opening angle model (ILamb et al.l 12006 



Donaghy 1 120061 ) . expect a positive correlation between E^ k and the bulk Lorentz factor. 
In the off-axis jet model (jYamazaki et al.l 120041 ) . no correlation is expected between 70 and 
-^pcak because the Doppler factor will change as a function of a viewing angle but not 70 . 
Another interesting theoretical model to discu ss is the case of a very high Lorentz factor 



( jMochkovitch et al.l 120031 ; iBarraud et al.l 120051 ). According to this model, XRFs can be 



produced in a condition with a very high 70 (so called "clean fireball"), while classical GRBs 
have a moderate 70 . In this case, we would expect a negative correlation between E^ k and 
the bulk Lorentz factor. Some additional effort, for example, to estimate -Speak from the 
Swift BAT data (Sakamoto e t al. in preparation) and to estimate a redshift from the Swift 



data (e.g., iGrupe et al.ll2006l ) is encouraged in order to discuss the correlation between 70 



and -Ep^ak an d the origin of XRFs. 

We can derive another constraint on 70 which is independent of the previous discussion. 
Our observational results suggest that the photons originating from an internal shock via the 
curvature effect (PLO component) and photons from an external shock (EXP component) 
arrive at the observer almost simultaneously. The observed time of the photon from an 



4 http://heasarc. gsfc.nasa.gov/docs/swift/archive/grb_table/ 
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internal shock, T inU can be expressed as, 

Tint ~ (Rint/2cj int )[l + (jintOf], 

where Ri n t is the radius where an internal shock emits, 7^ is the bulk Lorentz factor at Ri n t, 
and 9 is the jet opening half-angle. On the other hand, the observed time of an external 
shock emission, T ext , is expressed as, 



T ext ~ iW( 4 C7L), 



where R PX + is th e radius where an external shock emits and ^ exi is a bulk Lorentz factor at 
Rext ( ISarilll998l ). If we assume T int ~ T ext and also j int ~ ^ ext ~ 70, we have, 

2i^ nt [l + ( 7o #) 2 ] ~ i? 

-^peak can be written as a function of 6* and 70, 

3£k(*) ~ (27oH)[l + (lo0) 2 ] ~ [Elll k (9 = 0)]/(l + (7o#) 2 ), 

where Sp^^ = 0) * s S pea k observed by the on-axis observer. In this case, E^ k (9 = 0) 
corresponds to the observed S peak multiplied by (1+z). Therefore, the relationship between 
Rint, Rext and S peak is given by, 

R ext /2R int ~ [.Speak = 0)]/ E peak (9). 

Since the observed photon index of the XRT steep decay emission of our sample is ~ 2, which 
suggests that the observation of the spectrum is above S peak , it is reasonable to assume that 
the upper limit of S peak (#) is in the few keV range. Hence, the condition will be, 

Rext/^Rint > [-Speak = 0)] /S peak (#') 

where E ppa \r {9') is the upper limit of a few keV. If we use the angular spreading time (AT ang ) 
( lPiranlll999l ) as the time scale of an internal shock, then, 

Rint ~ 2c7 AT ang . 

From Eq. (j3J), the radius of an external shock can be expressed as 

R ext = (47rm p c 2 /3)- 1/3 (S /n /5M ) 1 / 3 7 - 2/3 . 

Thus, 70 can be written as, 

70 > (l/4) 3 / 8 (47rm p c 2 /3)- 1 / 8 (So/n /5M ) 1/8 (cAr an9 )- 3 / 8 [S P e ak (0 = 0)/S pcak (^)]" 3/8 . (7) 
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In the case of GRB 050803, AT ang can be derived as the duration of the last spike (the 
pulse at to+90 s), AT ang ~ 3 s, if one takes into account the time dilation effect. Although it 
is not possible to extract the information about the -E pea k of this pulse from the BAT data, 
it is reasonable to assume -E pea k > 100 keV since the photon index from a simple power-law 
fit to the BAT data is 1.2 ± 0.2 which is close to the low energy photon index of the typical 
GRB spectrum (-E pea k should be around or above the BAT upper energy range of 150 keV). 
Therefore, we use E pcak (6 = 0)/E pcak (6') > 100. The lower limit of 7 of GRB 050803 is 
estimated to be < 50 assuming a redshift of 2.4 (the mean redshift of the Swift long GRBs) 
and Eq/ujsm of 1 x 10 52 ergs cm 3 . Applying the same assumptions for the GRB spectral 
parameters (E peak {6 = 0) / E peak {6') > 100), we estimate a lower limit of < 60 for GRB 060418 
using the measured redshift of 1.489, AT ang of 2 s (the pulse at to+50 s), and Eq/uism of 
1 x 10 52 ergs cm 3 . The estimated lower limit of 70 for GRB 060418 is not contradicted by 
the value based on the REM observation. In summary, for the EXP component, one can 
provide a reasonable bulk Lorentz factor, 70, within the external shock scenario, which is 
consistent with other measurements. 



4-2.2. Internal shock scenario for EXP component 



As we discussed in §4.1, a simple interpretation based on a uniform jet model could 
contradict the prediction of the curvature effect. Here we discuss an intern al shock scenario 



for bo th EXP and PLO components based on an inhomogeneous jet model. lYamazaki et al. 



( 120061 ) investigated the GRB prompt emission 100-1000 seconds after the GRB trigger within 
the frame work of a multiple sub-shell model. According to their study, despite an angular 
inhomogeneity of the jet, the tail emission has a monotonic decay which resembles the 
XRT steep decay. In this context, if the jet has a core in which the emission energy is 
densely confined compared with the outer region, the PLO decay component arising from 
an on-beam sub-shell may be overlai d by th e off-beam core emission which causes the EXP 



decaying component. iTakami et all (120071 ) further extended their study, and in order to 



investigate the unknown jet structure, they proposed unique definitions of the decay index 
derived by unique definitions of the time-zero and of the fitting interval of the observed 
light curve. They found that the decay index in their definitions should have a wide scatter 
in the case of a power-law like structured jet. Here, we calculated t he decay index usin g 
our BAT and XRT composite light curve based on the definitions of ITakami et al.l (120071 ). 
Because of the diffi culty of using exactly the same definition of time-zero as proposed in 
Takami et al.l (120071 ) (T* in their paper), we define the end of the BAT emission as T*. We 
can fit our light curves using the proposed fitting interval for five GRBs in our sample (a ta u 
as a decay index, \ 2 an d d.o.f. of the fits are shown in the last two columns of table |3]). 
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Our results, based on a very small sample, show that the decay index ranges from 0.6 to 
3.0. However, both the size of the sample and the number of data points included in the 
light curve fit are very small because of using XRT data above 2 keV. There is an additional 
issue further reducing the sample; the appearance of a shallow decay in the XRT data in the 
fitting interval. Unfortunately, it is hard to conclude the structure of a jet with our limited 
sa mple. Once the GR B sample suitable for fitting a light curve with the unique definitions 
of iTakami et all (120071 ) can be increased, we may be able to draw a conclusion about the jet 
structure of GRBs using the XRT steep decay component. 



5. Summary 

In this paper, we presented a systematic study of the steep decay emission observed by 
the XRT. We constructed composite light curves in the 15-25 keV band extrapolating the 
XRT data (2-10 keV) up to the BAT energy range (15-25 keV). Based on the simultaneous 
fitting of the BAT and XRT data, we confirmed the existence of an EXP component for the 
majority of the bursts in the sample. We found that for the PLO component, the majority 
of the GRBs in our sample are inconsistent with the relationship of the curvature effect, 
a — r + 1, which is only valid in the case of the uniform jet. We also found that more than 
15% of the prompt fluence may be radiated below the BAT sensitivity limit for half of our 
sample. We argue that the EXP component could be the emission from the external shock 
which may indicate the deceleration of the initial shell by ISM during the prompt phase. We 
discuss the case of the prompt tail emission from the structured jet as an origin of the XRT 
steep decay but the sample is too small for a solid conclusion. 
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Table 1. Parameters of the light curve fits. Errors quoted at 68% confidence level. See text for details (section 2.3 

and 3). 
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3.5lo 8 x 10~ 9 
(3.4 ±0.3) x 10" 9 
2.8 X 10 _9t 


15.2 / 18 
118.2 / 27 
22.1 / 26 


060211A 


2006-02-11 
09:39:59.9 


XRT 
PL/EX 
PLE 


232-312 
0-913 
0-913 


ot 

72 
72 


2.1 ±0.9 
2.3 

2.29 ±0.02 


1.6 x 10~ 5 

1.7 x 10" s 
1.4 X 10~ 5t 


5.0 / 6 
43.8 / 29 


LZ < -39 

81 ±2 
54 ±4 


8.3+^ 2 x 10- 10 
(3.6 ±0.2) x 10" 9 
3.1 X 10- 9t 


5.2 / 6 
333.0 / 32 
39.9 / 31 


060306 


2006-03-06 
00:49:09.3 


XRT 
PL/EX 
PLE 


97-147 
35-256 
35-256 


ot 

40 
40 


3.7 ±0.8 
2.3 

2.33 ± 0.02 


2.2 x 10~ 3 
1.1 x 10~ 6 
1.0 X io~ 6 t 


3.1 / 2 
31.4 / 4 


3lt y 6 
17.8 ±0.7 
13 ± 1 


i.7t:to x 10- 9 

(3.1 ±0.7) x 10~ 8 
2.5 X 10~ 8 t 


3.3 / 2 
145.6 / 6 
25.4 / 5 


060418 


2006-04-18 
03:05:49. 2 b 


XRT 
PL/EX 
PLE 


178-400 
108-797 c 

108-797 6 


ot 

146 ± 2 
148.7 ±0.3 


2.6 ±0.1 
1.07 ±0.08 
0.89 ± 0.01 


3.1±t7 x 10" 4 
2.2l ( 1 ,;g x 10~ 8 
7.5 X 10~ 9 t 


27.1 / 28 
25.0 / 30 


101 ±6 
262 ± 23 
43 ±3 


(2.0 ±0.3) x 10- 9 
(3.0 ±0.5) x 10" 10 
1.5 X 10~ 8 t 


31.6 / 28 
123.0 / 31 
23.4 / 30 


060427 


2006-04-27 
11:43:01.0 


XRT 
PL/EX 
PLE 


148-198 
0-218 
0-218 


ot 

66 
114 


4.4 ± 1.3 
2.8 
6.9 


4.8 x 10" 1 
3.7 x 10" 6 
2.7 f 


3.7 / 3 
3.5 / 3 


41 +2U 
41 -10 
47 ± 2 

47 


5.lt 4 3 s x 10- 9 
(2.9 ± 0.4) X 10" 9 

2.8 x 10" 9t 


4.0 / 3 

5.1 / 6 

3.3 / 5 



060428B 2006-04-28 XRT 235-340 5 1 5.0 ± 0.7 74 7.9 / 7 55+i u 4.1+S'? x 10~ u 9.4 / 7 



Table 1 — Continued 



GRB 


to 


Data* 


Fitting 




Power-law 








Exponential 






UT 




[s] 


, pow 
t 


a 


-K"pow 




x7dof 


w 


^"cxp 


X^/dof 




08:54:15.2 


PL/EX 
PLE 


0-340 
0-340 


229 
236 + 2 


0.7 
0.40 ± 0.05 


3.5 x 10~ 
6.8 x 10" 


-10 

nt 


2.9 / 6 
2.7 / 7 


48 ±0.9 

39 ± 2 


8.9 X 10" 9 

1.0 x 10~ 8 t 


10.4 / 8 

2.7 / 7 


060923C 


2006-09-23 
13:33:10.8 


XRT 
PL/EX 
PLE 


205-529 
0-529 
0-529 


ot 

-3 ± 2 

-2.9 ±6 


2.7 ±0.3 
2.76 ± 0.02 

2.42 ± 0.03 


1.7 x 10 
2.0 X 10 

2.6 x 10" 


-4 

-4t 

-Bt 


3.5 / 5 
3.5 / 6 

3.2 / 5 


43 ±1 

oo+4 
oo -10 


4.0 x 10" 10 
(1.3 ±0.2) x 10~ 8 
9.5 x 10- 9 t 


5.4 / 5 
45.3 / 6 
3.2 / 5 



*XRT: fitting only the XRT data with a PLO model and an EXP model, PL/EX: joint BAT and XRT fit with a PLO model and an EXP model, 
PLE: joint BAT and XRT fit with a PLEXP model. 

tFixcd value. 

a All the BAT data points and the XRT data from to±224 s to t +350 s are used in the fit. 

b Although battblocks found a time interval which is 63 second before to, we concluded that this interval is due to the contamination of Sco X-l in 
the BAT field of view based on the BAT image analysis. 

c The fit to the XRT data is from t +108 s to to+148 s, and from t +797 s to t +400 s. The BAT data points from t +146 s to t +156 s are also 
included in the fit. 
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Table 2. XRT spectral parameters based on a joint fit to WT and PC data using data 
above 2 keV. The error in the photon index is quoted at the 90% confidence level. 



GRB 


WT Fitting Range 
[s] 


PC Fitting Range 
[s] 


FxRT 


x 2 /dof 


050803 


100.1- 


-184.1 


185.5- 


-12976 


1.9 ±0.2 


33.0 / 37 


050814 


166.6- 


-384.9 


386.5- 


-14133 


2.1 ±0.2 


33.5 / 35 


050915B 


158.4- 


-313.3 


229.9- 


-13791 


1.8 ±0.3 


9.3 / 11 


051109A 


131.4- 


-214.3 


3444- 


17540 


2.2 ±0.2 


19.0 / 31 


060109 


109.8- 


-200.2 


201.6- 


-12559 


2.0 ±0.2 


23.7 / 21 


060111B 


89.2- 


426.7 


155.4- 


-12655 


2.2 ±0.3 


10.0 / 14 


060202 


175.1 


-1025 


1028- 


13027 


i qq+0-05 


364.6 / 337 


060211A 


137.1- 


-330.9 


332.5- 


-13054 


1.8 ±0.1 


68.6 / 66 


060306 


96.9- 


174.5 


175.8- 


-12874 


2.40 ±0.03 


9.8 / 15 


060418 


166.1- 


-697.7 


699.9- 


-12561 


2 00 +0 - 07 


211.7 / 209 


060427 


148.2- 


-236.5 


237.9- 


12691. 


2.0 ±0.4 


7.7 / 10 


060428B 


235.6- 


-440.8 


442.7- 


-12773 


2.6 ±0.3 


16.6 / 17 


060923C 


204.6- 


-267.4 


268.7- 


-12796 


2.1 ±0.4 


3.0 / 7 



Table 3. BAT prompt emission properties, and the estimated Bulk Lorentz factor, 70, and the radius of the e xtern al 
shock, R based on the external shock model. The XRT decay index based on the definition of iTakami et al.l (120071 ). 
a t aii, and x 2 °f t ne fit are also shown. Errors quoted the 68% confidence level. 



GRB 


T90 

[si 


Model* 


TflAT 


-^pcak 

fkeVl 

L v J 






S3/S2 C 


z 


T±/2 


t 

To 


R 
Fcrnl 


a tail 


X 2 /d-o.f. 


050803 


87.9 


PL 


1.39 ±0.07 




22 ±1 


8.1 ±0.7 


1.5 ± 0.1 




18 


270 


2.8 x 10 16 






050814 


140.6 


PL 


1.8 ±0.1 




19 ±1 


6.2 ±1.2 


1.2 ± 0.1 


5.3 1 


45 


240 


3.0 x 10 16 






050915B 


40.9 


CPL 


1.4 ±0.2 


Mi; 


34 ±1 


17 ±1 


1.05 ± 0.04 




24 


242 


3.0 x 10 16 






051109A 


37.2 


PL 


1.5 ±0.1 




22 ± 2 


29 ± 3 


1.4 ±0.2 


2.346 2 


15 


287 


2.7 x 10 16 


0.59 ± 0.05 


9.7 / 14 


060109 


115.4 


PL 


1.9 ±0.1 




6.6 ±0.6 


3.4±0.1 


1.0 ±0.2 




23 


246 


3.0 x 10 16 






060111B 


58.8 


PL 


1.0 ±0.1 




16 ±1 


14 ±2 


2.0 ±0.1 




9 


350 


2.4 x 10 16 


0.9 ±0.1 


1.0 / 3 


060202 


198.9 


PL 


1.8 ±0.1 




22 ±1 


3.7 ±0.8 


1.2 ± 0.1 




98 


143 


4.3 x 10 16 






060211A 


126.3 


CPL 


0.9 ±0.3 


58+4 


16 ± 1 


3.3 ±0.1 


1.1 ± 0.1 




24 


242 


3.0 x 10 16 


1.3 ±0.2 


1.0 / 2 


060306 


61.2 


PL 


1.80 ±0.05 




21 ±1 


47 ±2 


1.1 ± 0.1 




46 


190 


3.5 x 10 16 


1.2 ± 0.2 


2.4 / 3 


060418 


103.1 


PL 


1.64 ±0.03 




80 ±1 


49 ± 2 


1.28 ± 0.03 


1.489 3 


21 


226 


3.1 x 10 16 






060427 


64 


PL 


1.9 ±0.2 




5.0 ±0.5 


1.7 ±0.7 


1.1 ± 0.2 




32 


217 


3.2 x 10 16 






060428B 


57.9 


PL 


2.6 ±0.1 




8.2 ±0.5 


3.4 ±0.6 


0.7 ±0.1 




27 


232 


3.1 x 10 16 


3.0 ±0.3 


4.7 / 8 


060923C 


75.8 


PL 


2.3 ±0.1 




16 ±1 


5.0 ±1.5 


0.8 ± 0.1 















*PL: power-law (dN/dE ~ E- r B^), CPL: cutoff power-law (dN/dE ~ E~ f bat exp(-(2 - T BAT )E / E veak )) 

t Calculated bulk Lorentz factor assuming Eo/n = 1 x 10 52 erg cm 3 . 

a Energy fluence in the 15—150 keV band [10 — 7 erg cm -2 ] 

b l-s peak energy flux in the 15—150 keV band [10 -8 erg cm -2 s _1 ] 

c Fluence ratio between S3(50-150 keV)/S2(25-50 keV) 

^Jakobsson et all ll2006r i 

2 Quimby et all 11200511 

^Dupree et al. I l l2006h 
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Fig. 1. — Schematic figures of the observed spectra with different -E pe ak energies (case 1: 
Epeak > 150 keV, case 2: 15 keV < -E pea k < 150 keV, case 3: 10 keV < E pcak < 15 keV, case 
4: 0.3 keV < -Speak < 10 keV, and case 5: E peak < 0.3 keV) demonstrating issues with the 
different extrapolations (left: extrapolating the BAT data down to the XRT energy band 
(BAT-to-XRT extrapolation); right: extrapolating the XRT data up to the BAT energy 
range (XRT-to-BAT extrapolation)). The light blue hatched regions show the extrapolated 
energy band. The extrapolated spectra shown with red dotted lines indicate issue with the 
extrapolation. There should be no issue with the extrapolation of the spectra shown in case 
1 and 5 (blue dotted lines). See text for details (section 1). 
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Fig. 2 



- The BAT and XRT composite light curve. See text for details (section 3). 
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Fig. 3. — continued. 
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Fig. 4. — continued. 
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Fig. 5.- 



continued. 
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Fig. 6. — continued. 
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Fig. 7. — The relationship between the power-law decay index, a, of the best-fit light curve 
with a PLO component and a photon index, Txrt- The dashed line is the expected rela- 
tionship from the curvature effect (a = 1 + Txrt)- 
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Fig. 8. — The relationship between the decay constant, w, of the exponential model and the 
BAT prompt emission properties. Top left: w vs. BAT T 90 , Bottom left: w vs. the ffuence 
ratio between the 50-100 keV and 25-50 keV band (the dotted and dash-dotted lines are the 
calculation assuming E pcak = 30 keV and E pcak = 100 keV, respectively, with a low and a 
high energy photon index of 1 and 2.5 in the Band function), Top right: w vs. the fruence 
in the 15-150 keV band, Bottom right: w vs. the 1-s peak flux in the 15-150 keV band. 
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Fig. 9. — A histogram of the ratio, in percent, between the fluence accumulated from the 
end of the emission as detected by the BAT to 1000 s after to (S (Tail: 15-25 keV)) and the 
fluence recorded by the BAT (S(BAT:15-25 keV)). 
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Fig. 10.— The BAT and XRT composite light curves of GRB 050803 (top) and GRB 050814 
(bottom) overlaid with the best-fit PLO model for different cho ices of a and tn° w ■ The labels 
in the parentheses are (tQ OW , a). The model which represents iLiang et al.l (120061 ) is shown 
with a blue dashed line. 



